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Stability of a Bose-Einstein condensate in a driven optical lattice: Crossover 
between weak and tight transverse confinement 
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We explore the effect of transverse confinement on the stability of a Bose-Einstein condensate 
(BEG) loaded in a shaken one-dimensional or two-dimensional square lattice. We calculate the 
decay rate from two-particle collisions. We predict that if the transverse confinement exceeds a 
critical value, then, for appropriate shaking frequencies, the condensate is stable against scattering 
into transverse directions. We explore the confinement dependence of the loss rate, explaining the 
rich structure in terms of resonances. 

PACS numbers: 67.85.Hj, 34.50.-s,03.75.-b 


I. INTRODUCTION 

Intense effort has been directed at using periodic 
driving to change the properties of cold atom systems 
Qi- For example, several experimental groups have 
been able to control effective hopping matrix elements 
by lattice shaking (MU- The Sengstock group explored 
the physics of classical frustrated magnets [Ifl, dll ■ The 
Ketterle group saw evidence for Hofstader butterfly 
physics [ij, The Esslinger group realized the 

topological Haldane model [I^ . The list of experiments 
and theoretical proposals is large [l5l - l53| . There are, 
however both conceptual and practical issues with using 
periodic driving to control a system. A driven system 
has no “ground state” nor a well-defined thermody¬ 
namic temperature. Furthermore, nearly all successful 
examples of this technique study non-interacting or 
very weakly interacting particles, and one almost 
always sees strong heating effects when moderate or 
strong interactions are introduced. Here we model 
some simple examples where these fundamental and 
practical issues are transparent. We make a series of 
predictions which are readily verifiable us ing techniques 
demonstrated in recent experiments |l5l . H^ . and 
which will enable the experimental study of interacting 
Floquet systems with cold atoms. 

In earlier works we began studying these questions by 
modeling two geometries. First, we considered a ID gas 
of atoms trapped in a shaken ID lattice [^. There we 
found large parameter regions where a Bose-Einstein 
condensate (BEC) is stable against 2-body collisions. 
Second, we considered a 3D gas of atoms trapped in a 
shaken ID lattice, making an array of “pancakes” 

We found that two-body collisions allowed energy to be 
taken from the shaking and transferred to transverse 
motion. The heating rates were consistent with those 
observed in experiment. In this setting, there is no 
steady state: the energy increases monotonically with 
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time. The natural question is how these limits are 
connected. A 3D gas with harmonic transverse con- 
hnement should interpolate between these behaviors. 
Here, we calculate heating rates in this crossover. 

We hnd a rich structure. First, there is a crit¬ 
ical strength of the transverse confinement beyond 
which two-body collisions are unable to deplete the 
condensate. Second, as a function of the transverse con- 
hnement, the dimensionless loss rate is non-monotonic, 
displaying drops and jumps characteristic of resonances. 
We explain this behavior in terms of the opening and 
closing of transverse decay channels. Our results will 
be crucial to the next generation of experiments. For 
example, one will be unable to observe a Floquet 
fractional quantum Hall effect without tuning to 
parameters where losses are negligible. While other 
authors have conducted related studies of the stability 
of driven systems the question of collisional 

loss into transverse channels is relatively unexplored. 
While we focus on a particular model, the loss into 
transverse modes is quite generic in cold atoms (see for 
example [HI). 

Driven systems are also of interest outside of cold 
atoms. For example, Photonics experiments have seen 
analogs of topological insulator physics in geometries 
which mimic electrons in a honeycomb lattice subjected 
to appropriate driving (^ . There are no interactions 
here. Translating these ideas to an electronic setting 
will require understanding the sort of loss processes 
which we study. Given the different structure of the in¬ 
teractions, our results cannot be applied directly in the 
solid state setting, but our logical framework is valuable. 

Figure 1(a) depicts a 1-D lattice with weak transverse 
confinement yielding an array of pancake traps. We 
consider driving the system by moving the lattice sites 
back-and-forth in the lattice direction. Figure 1(b) 
illustrate the tight confinement limit. Figure 1(c) illus¬ 
trates a 2D lattice in the weak confinement limit, where 
one has an array of cigar shaped traps. We consider 
square arrays, with the shaking oriented 45° from a 
lattice direction. These geometries are motivated by 
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the experiments performed at Chicago [mil- Using a 
kinetic model, we predict the scattering rate of bosons 
from a BEC as a function of the transverse confinement. 
Bilitewski and Cooper have performed a related study 
of the population dynamics in the Floquet realization 
of the Harper-Hofstader model [fi^. In places where 
our studies overlap our results agree. 

In section II, we introduce our model for analyzing 
the shaken lattice experiments. We also discuss the 
general formalism for obtaining the Floquet band 
structure. In section III, we use Fermi’s golden rule 
to predict the scattering rate for bosons out of the 
BEC and obtain the stability phase diagram for a 
BEC loaded in a one-dimensional shaken lattice and 
in a shaken square lattice. Finally, we conclude with 
directions for future experiments. 


II. MODEL 

Energy is not conserved in a periodically driven 
system since the Hamiltonian, H{t) is time dependent 
and work is being done on the system. However, one 
can map the problem onto a static one by viewing the 
system stroboscopically i.e. at times 0, T, 2T... nT (cf. 
ref. i)- At these times, the evolution operator is given 
by: 


U{nT) = C/(T)” =rexp dt H{t)/h 

where T stands for the time-ordering operator. This 
structure allows us to define a time-independent effec¬ 
tive Hamiltonian, Hgg via 



U(T) = 


Texp i J dt H{t)/h 


( 2 ) 


This defines the effective Hamiltonian, i/eff and the 
eigenvalues of the effective Hamiltonian, the quasi¬ 
energies. It is evident from Eq.Q that the effective 
Hamiltonian is not unique. If e is a quasi-energy, 
then e -I- mhio is also a quasi-energy where m is an 
integer and uj = 2ttIT. Clearly, the system has no 
unique thermodynamic ground state and the eventual 
occupations of various modes can only be determined 
by detailed modeling of the kinetics. 


For our model, it suffices to use the rotating wave 
approximations (RWA) and we do not need the full 
Floquet formalism. While this simplifies the mathe¬ 
matics, the conceptual issues are unchanged. 



FIG. 1. (Color Online) Schematic of shaken optical lattices: 
(a)ID lattice with weak transverse confinement; (b) ID Lat¬ 
tice with tight transverse confinement; (c) 2D lattice with 
weak transverse confinement. Ellipsoids represent edges of 
cloud in each well of the optical lattice sites and arrows il¬ 
lustrate motion of trap. A typical spacing between lattice 
sites is 532 nm (half the laser wavelength \l = 1064 nm) 
and a typical shaking amplitude is 15 nm. 


A. One-Dimensional Shaken Lattice 


The starting point of our modeling is the set-up 
in [l^ where a BEC of ^^^Cs atoms is trapped in 
a one-dimensional shaken optical lattice (with weak 
transverse confinement). When the shaking amplitude 
exceeds a certain critical value, the BEC undergoes a 
phase transition to a Z 2 superfluid (where condensation 
occurs at finite momentum k = ±fco yf 0). A schematic 
of the dispersion is shown in Fig. [21 For modeling this 
physics, it is sufficient to consider the first two Bloch 
bands and ignore the remaining bands (see supplement 

of [13). 


In the frame of the moving lattice, the Hamiltonian 
for the driven system is given hy H = Ho{t) -|- Hint, 
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FIG. 2. (Color Online) Schematic showing first (top) and 
second (bottom) Floquet quasi-energy bands of an optical 
lattice: e is the single-particle energy (arbitrary units used 
for schematic), k is the quasi-momentum and a is the lattice 
spacing. Since Floquet energies are only defined modulo 
the shaking quanta hw, the energy of the second band has 
been shifted down by hio so that it lies below the first band. 
Alternatively, this shift can be interpreted as working in a 
dressed basis, where the energy includes a contribution from 
the phonons. The mixing between the bands depends on 
the shaking amplitude. Dashed curves correspond to weak 
shaking, where the first band has its minimum at A: = 0. 
Solid curves correspond to strong shaking, where there are 
two minima at A: = ±A:o 7 ^ 0. 


where [I|, 


written as 

n,k 

+ Fo cos{wt) (xa^b^ + (4) 

Here, f-nki^nk) i® dispersion of the first (second) 
band, a"( 6 ") is the annihilation operator for particles 
in the first (second) band with the harmonic oscillator 
level being n and x is dipole matrix element between 
the first and the second band. As described in Ap¬ 
pendix A, Cnfc is generally time-dependent. However, 
when FQa/{hiu) <C 1 , Cnfc can be taken to be time- 
independent. 

We make the transformation 6 k exp(—ia;t)5k and 
discard far off-resonant terms (making the rotating wave 
approximation) to simplify the single-particle Hamilto¬ 
nian : 


^(sp) _ ,(1) nt n , (2),nt.n 

^RWA — 2 ^ ^nfc®fc + ^nfc'^fe "fe 
i.k 


XF (of fofc + 


(5) 


Here = ei^^ + K+% + 1)^^: = 4^^ + K + ny-b 

l)hfl — hio. We diagonalize this quadratic form writing 


rr(sp) _ 

-"rwa ~ 



( 6 ) 


For a particular value of n = {nx,ny}, the dressed 
dispersions and are shown as solid lines in 

Fig. [21 The bare dispersions and are shown as 
dashed lines. Reference illustrates how the Floquet 
dispersion can be calculated when the RWA is not 
applicable. 


-I- J 'I''l'(r) ( 2 F 0 cos(u;t))'l'(r) 

-I- J(fr ^''l'(r) 

Hint = 1 J(fr «'^(r)T'1'(r)^'(r)^'(r). (3) 

The atomic mass is m, the wavelength of the laser 
forming the optical lattice is Xl, the force from the 
periodic shaking is Fq cos{ojt) and g « is the 

interaction strength, being the scattering length. 
The transverse trap frequency is H. 


-h^ 

2m 


+ mVL^{x'^ -b y^) I fE'(r), 


B. Shaken Square Lattice 


We can easily extend the analysis of the previous 
section to the case of the square lattice. Since the 
shaken square lattice is separable and equivalent to two 
shaken one-dimensional lattices, one can write down 
the single-particle part of the Hamiltonian, H^o in the 
frame of the optical lattice as : 


H 2 D = J d^r^{ry{Hijj{z) + iJiD(y))4'(r) 


\ 2 m dx'^ 


+ mVt^x 


'I'(r)^f['(r) 


(7) 


As is detailed in Appendix A, the single particle part 
of the Hamiltonian describing the system iJgp can be 


+ 
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FIG. 3. (Color Online) Schematic showing the dispersion of 
the first Floquet band of a shaken square lattice beyond a 
critical amplitude. Color represents energy in units of the re¬ 
coil energy, Er (see scale). We see that the superfluid order 
parameter develops a D 4 symmetry in momentum space. 



where, 

(P Tr ■ 2 f 27r2:\ . 

= + (^_j+.FoCOs(a.t) 

Hio{y) = -b Vbsin^ (^) +y^ocos{ujt) 

( 8 ) 

Performing the same manipulations as in the last 
section, we end up with the following single particle 
Hamiltonian : 


rTRWA_ v- -(1.1) -nt -n ,^(1.2) T" 

.iJoP) - / 7,. I, ^h. U 1- Uu U Uu L> 


2D ^ ^ ^n^hx ^ky^kx ,ky^kx ,hy ^n^hx^ky^hx,ky'~^hx,ky 

n,ky,kz 


( 2 , 1 ) nt -n , p( 2 . 2 ) D" 

' ^n,kx,ky^kx,ky^kx,ky ' ^n,kx,ky kx,ky^kx,ky 


(9) 


where -b 7^^^ 

wneie — tn=0,kz ^ ^n=0,ky 


-b (n -b i/2)hn. 


Due to the separability of the square lattice, instead 
of the Z 2 reflection symmetry, the ground band devel¬ 
ops a D 4 symmetry for shaking beyond a critical force. 
We show this schematically in Fig. [3l 


III. STABILITY ANALYSIS 

In this section, we use a kinetic approach to inves¬ 
tigate the stability of a Floquet BEC as a function of 
transverse confinement. One would expect from our 
results in that the Floquet BEC would be stable 
if the transverse confinement exceeds a critical value. 
Eor both the ID shaken lattice and the shaken square 
lattice, we find this critical transverse confinement 
strength. Eor tighter potentials, the condensate is truly 
stable against energy-momentum conserving two-body 
collisions. We also identify several distinct signatures 


of interaction-driven scattering. 

Within our rotating wave approximation we can 
simply apply Fermi’s golden rule. The rate of scattering 
of two atoms out of the BEC is then given by: 

^ = f E - e.) (10) 

/ 

where 

(a°^)^ 

iV’/) = io) 

( 11 ) 

where, 'Ffc is a shorthand for representing 
{ofc, 6fc, Cfc, dfe},the state \tj}i) denotes the BEC where 
the bosons have condensed at momentum /cq. while 
I'l/'/) denotes a state where two bosons have scattered 
out of the condensate to momenta ko + k and kg — k 
respectively. The energies of the final states are e/ and 
ei respectively. If we did not use the Rotating Wave 
Approximation, a more complicated expression is nec¬ 
essary [ 2 ^. Using Eq. (nni we investigate the stability 
of a Floquet BEC. All our calculations are done for 
the experimental parameters of Ref. [l^ a lattice depth 
of 7Eh, where the recoil energy. Eh = h?/{ 2 m\\) 
where \l = 1064 nm and m = 133 amu. Eor these 
units, the zero-momentum bandgap for the ID optical 
lattice is is 4.96 Eh and the lattice is shaken at the 
blue detuned frequency of 5.5 Eh- It is reasonable 
to assume that loss is exponential. If not, Ea. dTUl) 
only describes the short-time behavior. At finite 
temperature, there are also heating processes involving 
one condensed atom and non-condensed atoms, or two 
non-condensed atoms. At typical BEC temperatures, 
these are negligible. 


A. One-Dimensional Shaken Lattice 

We first consider the case of a Floquet BEC loaded in 
a shaken ID lattice. Eor this case, the boson scattering 
rate in Eq. m can be expressed as : 

dN 27rg2A2 1 ^ 

It ~ T^LilE^ ^ ^ 

where F is the adimensional scattering rate, L is the 
linear system size and l± = ^Jh./{mfl). The detailed 
derivation and the expression for F are given in Ap¬ 
pendix B. F depends on the lattice depth, shaking 
frequency, shaking force and transverse confinement. It 
does not depend on the scattering length or the density. 
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FIG. 4. (Color Online) Adimensional scattering rate F as 
a function of the forcing amplitude, Fo in the limit of weak 
confinement into a ID lattice [Fig. 1(a)]. Blue, Dotted : 
h£l/E r = 0.04, Red, Dashed : hFl/ER = 0.08, Black, Solid: 
Analytic result from Ref. . 


Fig. |4] shows r vs Fq for weak transverse confinement 
{hQ/Efi = 0.04and0.08). For small Fq, F rises quadrat- 
ically and is roughly independent of O. For large Fq, a 
series of resonances are visible.The lifetime of the con¬ 
densate is given by : 


N mLl^a 
dN/dt ^ %halNT 


(13) 


Taking typical experimental parameters from the ex¬ 
periment in ref. jl5|. m = 133 amu, L = 30000 nm, 
Zj_ = 1000 nm as = 1.5 nm, N = 30,000 and F = 0.01, 
we get r ^ 1 s. 


To better understand the structure of resonances 
in Fig. m we plot F on a log scale as a function of 
the transverse confinement frequency, in Fig. [5l 
The set of vertical lines are given by the formula 
= Ei 2 l{ 2 n) and = £' 22 /( 2 n), where 

Ai 2 — e-TT + e_,r ~ 4eQ and E 22 — e-n- + e_ 7 r ~ • 

these are related to the bandwidths of the two bands. 
These energy values, E 12 and E 22 correspond to the 
maximum longitudinal energy transfer in two different 
scattering channels and the resonance structure in Fig. 
[5] corresponds to the closing of scattering channels. 
The factor of 2n corresponds to the spacing of parity 
allowed states. This structure can be understood by 
considering the energy and momentum conserving 
scattering processes in Fig. ([6]). The density of states is 
large when the final state has |fca| = tt. This resonance 
structure leads to special parameters where the BEG 
would be particularly stable or unstable. These reso¬ 
nances are a useful fingerprint of the loss mechanism 
and can be used in an experiment to test our model 
of interaction-driven instability. The dashed line in 
FigIS] corresponds to a lifetime of r Ri 10s (using the 
parameters below Eg. (1131) ). There is a large window 
around hfl ~ l.lEji, where the lifetime exceeds 10 
s. The BEG is completely stable against collisions 


FIG. 5. (Golor Online) Logarithm of the adimensional scat¬ 
tering rate, F in a ID lattice [Fig. 1(a),(b)] as a function of 
the transverse trapping frequency, D for a fixed value of the 
forcing amplitude, Fo = Fc, where Fk is the amplitude where 
the dispersion of the ground band is quartic near k = 0. 
Red vertical lines denote resonances at corre¬ 

sponding to the closing of scattering channels (see text).The 
black dashed line shows the value of ln(r) for different val¬ 
ues of the transverse confinement for which the BEG lifetime 
is greater than 10 s (assuming the parameters quoted after 

Eq.Hai). 



FIG. 6. (Golor Online) Schematic illustrating conservation 
of energy and momentum in two-body collisions in a shaken 
pancake lattice. Black dot denotes condensate in first band 
at fc = 0. Solid lines show first and second with no trans¬ 
verse excitations. Arrows denote an energy and momentum 
conserving collision. The resonances in Fig. 4,5 correspond 
to the situation where the final states have |fca| = tt 


hil > 2.05Eii. In Fig. [71 we show how the stability 
boundary varies with drive frequency. In terms of the 
dispersions of the two bands, the critical confinement is 
given by : 


nn = £(.2) + e(_2) _ 2eW (14) 

For larger 17, energy and momentum can’t be conserved 
in 2-body collisions. 
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FIG. 7. (Color Online) Stability phase diagram for a BEC in 
a driven ID lattice for a fixed value of the forcing amplitude, 
Fo = Fc- Here lu—uiq is the detuning of the shaking frequency 
ui from the zero-momentum bandgap cuq. 
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FIG. 8. (Color Online) Adimensional scattering rate F as 
a function of the forcing amplitude, Fo in the limit of weak 
conhnement {Ml = 0.0875^) for a 2D lattice [Fig. 1(c)] 



FIG. 9. (Color Online) Logarithm of the adimensional scat¬ 
tering rate F in a 2D lattice as a function of the transverse 
confinement, 57 for a hxed value of the forcing amplitude, 
Fo = Fc. The black dashed line shows the value of ln(r) for 
different values of the transverse confinement for which the 
BEC lifetime is greater than 10 s (assuming the parameters 
quoted after Eg. (IT6ll 1. 


The lifetime of the condensate is given by : 

dN/dt ^halNT ^ ’ 

Now, taking typical experimental parameters from the 
experiment in ref. [l5l| . m = 133 amu, Ly = 30000 nm, 
Lz = 30000 nm, lo_ = 1000 nm as = 1.5 nm, N = 30,000 
and r = 0.4, we get t ~ 0.73 s. 


B. Two Dimensional Shaken Lattice 


In this section, we explore the stability of a Bose- 
Einstein condensate loaded in a two-dimensional optical 
lattice. The stability analysis is very similar to that of 
the shaken ID lattice. The scattering rate of bosons can 
be written down as : 


dN _ 2 tt 1 

dt hi LyLzl± 


(15) 


where F is the adimensional scattering rate and 
and Ly denotes the linear system size in the z and y 
directions. The detailed derivation and the expression 
for r are given in Appendix B. 


The scatter of points in Fig. [5] is related to the 
resonances. Again, these can be explored by fixing 
Fq to some value (here, F^) and then plotting the 
scattering rate, T as a function of the transverse 
confinement, 51 as shown in Fig. [HI The black dashed 
line again corresponds to a lifetime of 10 s. There 
are specific value of 51 at which the scattering rate 
drops significantly. These values of 51 are shown 
as vertical lines in Fig. IH] and correspond to hfl = 
{Ei 2 + E 22 )/ (2'n.), {Ei 2 + Ei 2 )I (2n), {E 22 + E 22 )/ (2n). 
As in the ID case, these frequencies correspond to the 
closing of scattering channels. There is also structure 
related to the van Hove singularities in the density 
of states, but for clarity, we do not mark them with 
vertical lines. Beyond a transverse confinement of 1.4 
Er, the BEC will almost always have a lifetime r > 10s. 


We show r for a relatively weak value of transverse 
confinement {HQ/Er = 0.08) in Fig. [51 We see that 
the adimensional scattering rate, T, is higher for the 
shaken two-dimensional square lattice when compared 
to the one-dimensional lattice. 


Due to the separability of the Hamiltonian, the 
critical transverse confinement for the 2D square lattice 
is exactly twice that of the one-dimensional lattice, so 
the stability phase diagram is readily inferred from FigjT] 
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IV. CONCLUSION 


where, 


In this paper we studied the effect of transverse 
confinement on the stability of a Floquet BEC for both 
a shaken ID lattice and a shaken 2D square lattice. 
We obtained scattering rates as well as the stability 
phase diagrams for both systems. The scattering rate 
shows a resonant structure and fine tuning parameters 
can drastically reduce the loss rate. This structure 
arises from the opening and closing of loss channels 
corresponding to the quantized transverse modes. It 
provides a fingerprint of the loss mechanism and could 
be a valuable tool for minimizing loss. We find a 
critical value of transverse confinement, beyond which 
there are no allowed 2-body scattering processes which 
can deplete the condensate. Well before this point 
however, the scattering rate drops to extremely small 
values, making the BEC stable for the time-scales of 
the experiment. 

The loss mechanism that we study has another dis¬ 
tinct signature - namely that energy is converted from 
the time-dependent potential into transverse motion of 
the atoms. This transverse motion can be directly 
probed in time-of-flight experiments. 
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Appendix A: Derivation of the ID Hamiltonian 


In a tight-binding prescription, the single-particle 
Hamiltonian describing the system in the frame co¬ 
moving with the lattice can be written as 


Xj ~ J zwl{z — Zj)w2{z — Zj) 

= j wliz-Zi) + 

t[f = J dz w;{z-z,) + 

F{t) = Fq cos(wt) (A3) 

Here, Wi is the Wannier function for the ith band. It 
should be noted that Xj is independent of j and so 
we can call it y. The operators aj and bj annihilate 
particles in the two bands. If necessary more bands can 
be included. 


Performing a basis rotation : |'0) —>■ Uc{t)\ip) where 


Uc(t) = exp ZjFo cos(wt)(a]aj -|- b]bj) j , 

(A4) 


we transform the Hamiltonian as: 


i7'(t) = UMt)U-^ - ihUcdtU-^ 

= E 


(A5) 


71 a; 


with 

'ij 

+ Y^Fcos{ujt) (xafbj + x*bfa]^^ 

3 

+ ^ h£l{nx + ny + l) 

n 

= J2J2^os{mka) (-Ji^^(t)af 

k TYl 

+ E ^0 cos(a;t) (yaf -k X*bl'^al) 

k 

+ Ml{nx + ny + l) + bfbl^ 

n 

(A6) 


Ho{t) = f d^rx_ Y + t?^b\hj + h.c^ 

+ E 

3 

+ -k j/^) + b^jbj^ (AI) 

(A2) 


where, 

J5(t) = exp(-*Fo^^(z, - zy)) 

= iy exp(-iFo^^4^a(i _ j))^ (A7) 

a = Al/ 2 is the lattice spacing and y = y* for a 
suitable choice of phase for ak and bk- We use n as a 
shorthand for denoting {nx,ny}. 






In the limit of Foa/{hu;) <C 1, Hence, we 

can write down the Hamiltonian as : 

^sp — 2^ Cjifc af. + o^, o^, 

n.k 

+ Fo cos{ujt) (xa^b^ + X^^fe^a") (AS) 

where, 

e^ = EE cos(mfca) + hFl{nx + + 1 ) 

k m 

e^ = EE cos{mka) + Ml{nx + riy + 1) 

k m 

(A9) 


Appendix B: Derivation of the scattering rate 
1. ID Lattice 


with 


r = 


L^lEna^ 

27r 


J dkrl‘^’’^^6{ef - e,) 


This is Eg. dT^ in the main text. 

2. 2D Square Lattice 

For the case of the 2D square lattice, the scattering 
rate in Eq. can be written down as 

^ = (B4) 

^ ' Ua,Tlh 

where 


For the case of the ID optical lattice, the scattering 
rate in Ea. dlOH can be written down as 


dN 2 tt ^^2 A 
dt h A 


E (HI) 


^ k 


and 




L-k^L+k^ko^'kol^z) 


LyLJ± 


where 


pna.rib 

^ k 


and 




Lll 


with na(n6) = {<,n^}({<,n^}), = 

7J„(a:) exp(—a;^/2), Hn{x) being the Hermite polyno¬ 
mial of order n. An important consequence of the form 
of lB2l is that = 0 unless and (as well as 

and n^) have the same parity. Finally, Eq. m can be 
simplified to write 


_ / dx(j)'^'^‘‘\x)(j)^'^'’\x)(j)^^\x)(j)^^\x) (B5) 


with l± = y/h/irn^ just as in the case of the ID shaken 
lattice. An important consequence of the form of IBSl is 
that is 0 unless ria and rif, have the same priority. 

Ean. (IB4l) simplifies to give : 


dN 27r 1 


dt h A LyLzl± Ena^ 


(B6) 


with 


- e.) 


dN 2 tt g^ fV2 I 


dt h A Llf Efia^ 


(B3) This is Eg. dTSl) in the main text. 
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